

rho_u = [0.1,0.8,0.95,0.99];
%rho_u = [0.1,0.99];

load ../toDynare
LL = 150;
X=1:1:LL;

% I renormlize to get the same NPV


for i=1:length(rho_u)
    name = ['fig_opti',num2str(rho_u(i)),'.mat'];
    load (name,'Bt_eps','ut_eps'); 
    if i==1
        rho_u1 = rho_u(i);
        eps = ut_eps(1);
    end
    coef = (1/ut_eps(i))*(1-rho_u(i)/(solution.R))/(1-rho_u1/(solution.R))
    %name = ['fig_opti',num2str(rho_u(i)),'.mat'];
    %load (name,'Bt_eps'); 

    plot(X,Bt_eps(1:LL)*coef,'LineWidth',2);




    hold on

end
legend(strcat('\rho_{\beta} =  ',num2str(rho_u(1))),strcat('\rho_{\beta} =  ',num2str(rho_u(2))),strcat('\rho_{\beta} =  ',num2str(rho_u(3))),strcat('\rho_{\beta} =  ',num2str(rho_u(4))),'Location','southeast');
hold off



legend boxoff 

